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Motivated by the recent interest in cosmologies arising from energy density modifications to 
the Friedmann equation, we analyse the scaling behaviour for a broad class of these cosmologies 
comprised of scalar fields and background barotropic fluid sources. In particular, we determine 
the corresponding scalar field potentials which lead to attractor scaling solutions in a wide class of 
braneworld and dark energy scenarios. We show how a number of recent proposals for modifying the 
Friedmann equation can be thought of as being dual to one another, and determine the conditions 
under which such dualities arise. 
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I. INTRODUCTION 

The Friedmann equation is one of the cornerstones of 
modern cosmology relating the expansion of the Uni- 
verse to the total energy density within it. It forms the 
starting point for almost all investigations in cosmology. 
However, over the past few years, possible corrections to 
the Friedmann equation have been derived or proposed 
in a number of different contexts, generally inspired by 
braneworld investigations. These corrections are often of 
a form that involves the total energy density p, and are 
such that they tend to play a role early in the history 
of the Universe, fading away as we enter the late-time, 
post-nucleosynthesis era (although that is not always the 
case as we shall see). Up to now, the different models 
have been presented in the literature without any at- 
tempt to relate them. In this paper, by introducing a 
generalised form for the correction, we will provide a for- 
malism which allows us to relate a large class of modified 
Friedmann cosmologies. Assuming the total energy den- 
sity to be comprised of a canonical scalar field <f> with po- 
tential V(4>), together with some form of barotropic fluid, 
we will demonstrate how the existence of scaling solutions 
determines the form of V((f>), and in doing so we will es- 
tablish a direct relation between the form of the potential 
and the functional form of the modification to the Fried- 
mann equation. Scaling (attractor) solutions in cosmol- 
ogy are very important because they allow one to under- 
stand the asymptotic behaviour of a particular cosmology 
and to determine whether such behaviour is stable or not. 
They have also been advocated by a number of authors as 
a way of establishing the behaviour of general scalar fields 
in a cosmological setting both in the context of conven- 
tional Friedmann cosmologies [JEHSSE0|l,H[l3 
and in particular classes of modified Friedmann cosmolo- 
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In section II we present the equations of motion aris- 
ing out of the modified Friedmann equations and in- 
troduce variables which allow the scaling solutions to 
be determined. The general conditions for scaling be- 



haviour are then established in Section III and we show 
that these can be written as a closed form relationship 
between the scalar field and the functional form of the 
modification to the Friedmann equation. In section IV 
we demonstrate the existence of duality symmetries be- 
tween different scaling solutions, and determine the con- 
ditions which must be satisfied in terms of the modifica- 
tion to the Friedmann equation for such duality proper- 
ties to be obtained. Section V then applies the results of 
the previous sections to a general class of models, eval- 
uating the scaling potentials (and their duals), as well 
as the explicit evolution of the scale factor in the scal- 
ing regimes. In particular, we apply the technique to a 
number of recently investigated cosmologies of relevance 
both to braneworld and dark energy scenarios, including 
the Randall-Sundrum |24|, Shtanov-Sahni [2j| and Car- 
dassian [2(| models. The Dvali-Gabadadze-Porrati [27| 
braneworld scenario is investigated in Section VI and we 
summarize our results in section VII. Throughout units 
are chosen such that % = c= 1 . 



II. EQUATIONS OF MOTION 

We consider spatially flat Friedmann-Robertson- 
Walker (FRW) cosmologies such that the dynamics is 
determined by an effective Friedmann equation of the 
form 



H 2 = 



8tt 
3mf 



pL 2 ( P ), 



(1) 



where H = a /a is the Hubble parameter, a is the scale 
factor, p is the total energy density of the universe, a dot 
denotes differentiation with respect to cosmic time and 
TO4 is the four-dimensional Planck mass. Modifications 
to standard relativistic cosmology are parametrized by 
the correction function L{p) and this is assumed to be 
positive-definite without loss of generality. 

We will investigate models where the universe is 
sourced by a self-interacting scalar field <f> with poten- 
tial V(4>) together with a barotropic fluid with equation 
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of state P 7 = (7 — l)p 7 , where 7 is the adiabatic index. 
The energy density and pressure of the scalar field are 
given by p^ = 4> 2 /2 + V and P$ = (j) 2 /2 — V, respec- 
tively. As in conventional cosmologies, we assume that 
the energy-momenta of these matter fields is covariantly 
conserved and this implies that 



p 7 = -37ifp 7 , 
4> = -ZH(j>-dV/d(j). 



(2) 
(3) 



Eqs. JIJ-JSl close the system that determines the cosmic 
dynamics. 

In standard cosmology the stability of scaling solutions 
is analyzed by introducing the variables 



X, 



CLW 



^CLW 




3m\ H 



(4) 



and rewriting the field equations as an autonomous sys- 
tem. Following [T2I IT^j , we define the new pair of vari- 
ables: 



X 
Y 



Vv 



that are related to those of Eq. Q by 
X Y 



X 



CLW 



Y< 



= eL, 



CLW 



(5) 



(6) 



where e = ±1 for expanding and contracting universes, 
respectively. In what follows, we consider expanding 
models unless otherwise stated. 

When expressed in terms of the new variables JSJ), the 
equations of motion Q JSJ can be written in the form: 



X' = -3X + eXJ^Y 2 + ^X[2X 2 + 7(1- X 2 - Y 2 ) 



y' = -e\^-XY+ i -Y[2X 2 + 1 {l-X 2 -Y 2 )], 



A' = -eV6A 2 (r-l)X 



where 



-3\[2X 2 + 7 (1 -X 2 - Y 2 )]p-[lnL(p)}, 
1 m4 dV/d(f> 



X 

r 



v 



V ' 

d 2 V/d<p 2 
(dV/dcpy ' 



(7) 

(8) 
(9) 

(10) 
(11) 



and a prime denotes differentiation with respect to the 
logarithm of the scale factor, N = In a. As in Ref. (lflj . 
Eqs. (fTTTft and (fTTT generalize the expressions introduced 
111 Refs. HSIlIli In particular, A is related to the 
parameter Aclw = — -j= dv ^ introduced in Ref. 
such that 



A 



Aclw 



1 

L' 



(12) 



We refer to A as the 'scaling parameter'. 

The system of equations Q- flU for the variables X, 
Y and A do not appear to be closed due to the presence 
of the term involving T in Eq. JHJ. However, it follows 
from Eqs. JSJl, (JTJJ and {TTJ that we have T = T (</>), 
(j) — (j)(p, Y) and p — p(X, 4>) and hence that V — T(X, Y). 
This implies that the equations are indeed closed. 

The definition of the total energy density implies that 
the variables JHJ satisfy the constraint equation 



X z 



Y 2 + ^ = l 
P 



(13) 



and, since the energy density of the barotropic fluid is 
semi-positive-definite, any cosmological model can be 
represented as a trajectory in the (X, F)-plane that 
is bounded within the unit circle, i.e., = p^/p — 
X 2 + Y 2 < 1. Furthermore, since Y > by definition, it 
is sufficient to consider the evolution in the upper half of 
the disc. 

Eqs. H7I9|I exhibit an important property. For the case 
where A is constant, they have an identical form to that 
of the plane-autonomous system of standard relativistic 
cosmology that is formulated in terms of the variables 
{Aclw j ^clw, Aclw}- This duality immediately implies 
that the system I|7I8|) admits an identical set of critical 
points to that of the standard scenario when these solu- 
tions are expressed in terms of the variables {A, Y, A}. 

A further consequence of such a duality is that the 
stability of each fixed point solution can be determined 
directly from the stability analysis of Ref. In to- 

tal, there are five critical solutions to Eqs. J7J) and © 
where the variables {A, Y, X} — {A c , Y c , X c } are con- 
stants. Three of these represent the unstable solutions 
(A c = 1, Y c = 0),(A C - -1, Y c = 0), (A c = 0, Y c = 0) 
for all values of A and 7. The value of A determines the 
nature of the other two points. For A 2 > 37, there exists 
an attractor solution 



Y r 



3(2-7)7 



2A 2 



(14) 



where the effective adiabatic index of the scalar field, 
defined by 



2X' 



70 



A 2 + Y 2 ' 



(15) 
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satisfies the condition j$ = 7. For this late-time attrac- 
tor solution, the relative contribution of the scalar field's 
energy density to the total energy density of the universe 
is constant, Vl^ c = X 2 + Y 2 = 37/A 2 , and consequently, 
the energy densities of the scalar field and fluid rcdshift 
at the same rate as the universe expands. 

The fifth critical point arises if A 2 < 6 and is given by 



X c = 



Yr 



_A_ 

v/6' 



1 



A 2 
~6~' 



(16) 



This corresponds to the case where the scalar field dom- 
inates the fluid (f^ c = 1) and has an effective adiabatic 
index 7^ = A 2 /3. The solution is stable if 7^ < 7, i.e., 
A 2 < 3 7 . 



III. GENERAL CONDITIONS FOR SCALING 
SOLUTIONS 

Since the expressions (01 and (|12|l relating the stan- 
dard and modified FRW cosmologies involve the correc- 
tion function L(p), the scalar field potential that gives 
rise to the fixed point attractor solutions (|14|l and i|16|) 
in a given generalized scenario will depend on the spe- 
cific form of this function. In particular, the potential 
will differ from the purely exponential form that leads 
to scaling solutions in the conventional FRW model. In 
this section we establish the correspondence between the 
modified Friedmann equation and the scaling potential. 

It can be shown by direct substitution that both sets 
of critical points 114(1 and Ijl6|l represent solutions to the 
field equations Q-© of the form X' = Y' = A' = if 
the relation 



T=l + p—[\nL(p)} 



(17) 



is satisfied. Since p = V/Y 2 for these solutions, Eq. fTTjl 
may be written in the form 



d 2 p/dcp 2 d 



(18) 



and multiplying Eq. (|18f) by (dp/ dip) / p then implies that 



d(\n(dp/d<p)) d(lnp) d(luL) 



dcj) 



0. 



(19) 



Eq. (|19fl may be integrated twice to yield a necessary 
and sufficient condition on the scalar field potential if 
the solution is to represent a scaling solution for a given 
choice of correction function L(p). We find the important 
result: 



dp 



8ttA 



ITi4 



(20) 



where one of the integration constants has been set to 
zero without loss of generality by performing a linear 
shift in the value of the scalar field and the constant of 
proportionality on the right-hand side follows by requir- 
ing consistency with Eq. H10|) . 

It is also of interest to determine the evolution of the 
scale factor for a given class of scaling solutions. Since 
X c is a non-zero constant for these solutions, Eq. (J5J 
implies that the scalar field is a monotonically varying 
function of proper time (0^0). It is natural, therefore, 
to view the value of the field as the dynamical variable of 
the system and to express all time-dependent parameters 
in terms of this variable. 

In general, the scalar field Eq. J21 can be expressed in 
the form 



or, equivalently, as 



-3H<P 2 



= -3H<j). 



(21) 



(22) 



It then follows from the definition of the Hubble param- 
eter that 



3H 2 = lda d P<t> 



a d<j) d<p 



(23) 



and substituting Eq. (|23|l into Eq. Q implies that the 
Friedmann equation can be expressed in the form 



da dp 



Sn 



dip dip fl^cm 2 
Introducing a new variable 



apL 2 {p). 



b(<p) = exp 



4>c 



dp 



pL 2 (p) 



simplifies Eq. to 

da db 



Sn 



ah 



(24) 



(25) 



(26) 



and the scale factor is then determined up to a single 
quadrature: 



a((p) — exp 



8tt 



d<pb((P) 



db 



(27) 



Thus, for a given cosmological scenario characterized 
by a correction function L(p), the potential (and equiva- 
lently the total energy density) yielding the scaling solu- 
tion is determined by integrating Eq. (|20|l . Integration of 
Eq. Ij25(l then yields the dependence of b(<p>) on the scalar 
field and the evolution of the scale factor follows after 
integration of Eq. I|27|l . Finally, the time-dependence of 
the scale factor can in principle be deduced by integrating 
Eq. (E2J|, 



t = 



dcPL^p 1 ' 2 ^) 



dp 



(28) 
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and inverting the result. 

In the following section, we employ the above formal- 
ism to establish a link between different classes of scaling 
solutions that arise for various choices of the modification 
to the Friedmann equation. 



IV. DUALITY BETWEEN SCALING 
SOLUTIONS 

A duality between different scaling solutions can be 
established by noting that Eq. l|26[l is invariant under 
the simultaneous interchange 



a(<t>) - V(4>), b{4>) - a 1 /? (0), 



(29) 



where p is an arbitrary constant. This symmetry im- 
plies that a given scaling solution may be employed as 
a seed to generate a new scaling cosmology for a differ- 
ent Friedmann equation and associated scalar field po- 
tential. To be specific, let us consider the scaling solu- 
tion parametrized by the functions {a+ (</>), b + (0), p+((j>)} 
that arises for a specific choice of correction function 
L + (p). We now denote the 'dual' scaling solution as 
{a_(0), b- (</>), p_ (4>)} and assume an ansatz of the form 

&_(*) = M0] 1/p - (30) 

The new scale factor is then determined from Eq. I|27|) : 



a_ (<j>) = exp 



8 ^P J d(f>a+W 



(31) 



However, since the function a+{<p) is itself a solution to 
the Friedmann equation Q26J1. Eq. (|3I(I simplifies after 
integration to 



o_(0) = [6+(0)]" 



(32) 



modulo an arbitrary (constant) prefactor. 

We may now determine the condition that the dual 
correction function L-{p) must satisfy for the solution 
to also represent a scaling solution that satisfies Eq. 
(|20|l . If we assume a priori that the two solutions a± (</>) 
represent scaling solutions characterized by A±, respec- 
tively, Eq. lj2T?|l implies that 



1 



dp- 



X+p+L+ dcf) \-p-L- dcf) 



(33) 



It then follows, after substitution of Eq. (|24|) into the 
right-hand side of Eq. (|33J) . that 



f dp^ 



8ttAi 



p+L + d<f) f^m^A- (da-/d(f>) 
and Eq. (|26[) then implies that 
1 dp + 



A+ L_ db- 



(34) 



(35) 



On the other hand, substituting the ansatz l)30f) into 
Eq. H35|) and employing Eq. (|24|l for the positive-branch 
solution yields the condition 



( ^f± 



8tt X a 



2 X 
4 A_ 



-L_L 4 



(36) 



Consistency with Eq. (|20|) therefore implies that a nec- 
essary and sufficient condition for the dual cosmology 
{a_(0), L_(0)} to represent a scaling solution is that 
the correction functions arising in the respective Fried- 
mann equations must be proportional to the inverse of 
each other when both are expressed as functions of the 
scalar field: 



(37) 



It is interesting that the standard relativistic cosmol- 
ogy (L = 1) represents the self-dual model when 

P n%x + X- = -l. 

Finally, we find after substituting Eq. (|3T)l into Eq. 
(|3*3^l and employing Eq. l(25*|l that the energy density of 
the dual scaling solution is given by 



1 dp- = _ n A 2 1 db+ 
p_ d(f> c ~ b + d(f> 



(38) 



The dual potential then follows immediately from Eqs. 
© and 123): 



V-{(j>) = F c 2 exp 



- P x 2 _n 



4>c 



dp^ 



(39) 



In the following section we employ the techniques de- 
veloped above to determine scaling solutions (and their 
duals) in a number of different cosmological settings. 



V. UNIFICATION OF MODIFIED FRIEDMANN 
COSMOLOGIES 

A. A Generalized Class of Scaling Cosmologies 

A wide class of scenarios that have been considered 
recently predict deviations from the standard cosmology 
of the form 



L(p) = ^1 + Ap-, 



(40) 



where v is an arbitrary dimensionlcss constant and A is 
an arbitrary constant with dimension m~ 4,y . In this case, 
the form of the potential leading to scaling (fixed point) 
solutions is determined by integrating Eq. (|20|) . It is 
found that 



V+(<f>) = ifA-^cosech 



if A > and 

V-(<f>) =Y?\A\- 1 /»sech 2/l/ 




Xv V8^f 

2 7Tt4 



(41) 



(42) 
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if A < 0. 

Given the form of the scalar potential (|41|l . the para- 
metric solution for the case A > is determined by inte- 
grating Eqs. fZty and (|2"7ll . respectively: 



p+{4>) = A- 1 / !/ cosech 2/l/ 



6+0) = A- f W l/ sech 2f W l/ 



a+(<£) = sinli 2 ^ 2 ^^ 




(43) 



The time dependence of the solution follows by substi- 
tuting Eq. (|4^|) into Eq. (|4T]|) to yield the Friedmann 
correction function: 



L + (<j)) — cotanh — 




and then evaluating the integrand in Eq. I|28|) : 



t = 



y/ZA 1 /" 



xn 



4>c 



.sinh 1 ^ - 




(44) 



(45) 



The integral l|45l) can be performed analytically for var- 
ious choices of v, whereas the late-time behaviour can 
be analysed for arbitrary v. In particular, we find from 
Eq. I|45|) that the late-time limit corresponds to large (f>, 
and it therefore follows from Eq. (|44|l that L — > 1 as 
t — > oo. This in turn is the limit corresponding to the 
case of an exponential potential. 

The corresponding scaling solution for A < driven 
by the potential l|42l) is deduced by applying the duality 
transformation (|29|l to the solution l|43[) for a particular 
value of the constant p, where the scaling parameters are 
chosen to be equal, A + = A_ = A. For the case where 
p = — l/(A 2 f2 2 c ), the duality transformation l|37|l implies 
that the dual correction function is given by 



L_ (4>) = tanh - 




(46) 



whilst integrating Eq. 1|39|) with the form for p+{<fi) given 
in Eq. I|43l) implies that the dual potential has precisely 
the form of Eq. I|42|) . We may conclude, therefore, that 
the dual correction function is given by Eq. (|40|l with 
A < 0. In this sense, a model with A > and a specific 
value of v is twinned with the model where the value of 
v is the same but the sign of A is changed. In general, 
the dual scale factor is deduced from Eqs. (|32l) and (|43ll : 



o_fo) =cosh 2 /( A2 °^ 




2 TO4 

and the time-dependence follows from Eq. (|28|l : 



t 



xn 



2 m<4 



(47) 



(48) 



It is of course trivial to show that for this choice of 
p, the standard cosmology solution corresponding to the 
case L = 1 reproduces the well known exponential po- 
tential for the scalar field [28j . 

In the following subsections, we consider some of the 
specific models that belong to the class of corrections 
given by Eq. |gU|) . 



B. Randall- Sundrum Type II braneworld 
cosmology 

The case A = 1/2<t and v = 1 corresponds to the 
Randall-Sundrum type II (R-S II) braneworld scenario 
pil l36l l3ll l38| . where a co-dimension one brane with 
positive tension a is embedded in five-dimensional Anti- 
de Sitter (AdSs) space: 



L{ P ) 



1 



P_ 

2(7 ' 



(49) 



For this case, the scaling potential yielding the fixed 
point solution is given by 



V(</>) = 2crY c 2 cosech 2 




(50) 



and the time-dependences of the scalar field and scale 
factor are deduced by evaluating the integral (|45|l for 
v = 1 and substituting the result into Eq. (|43l) : 




a(t) 



3Amj 



X fl^ct 



(51) 
(52) 



These results for the scaling solutions associated with 
the R-S II model confirm those previously obtained us- 
ing a different method in Ref. [33. An important fea- 
ture of our approach is that it shows the solution is a 
fixed point attractor solution. Another important feature 
that emerges is that in the limit where the quadratic en- 
ergy density term dominates, i.e., when \/8TTX(j) <C 2rn4, 
the potential asymptotes to V oc (f>~ 2 , consistent 
with earlier analyses |29l l30| . Similarly, once the energy 
density has decreased so that L ~ 1, the value of the 
scalar field becomes large, and V ~ exp[— X(y/8n /rni}4>], 
in agreement with general relativistic results [jj Q . 

In Figs. 1 1121 we confirm numerically how the above 
potential leads to the expected attractor solutions for a 
model with a barotropic fluid of radiation (7 = 4/3). 



C. Shtanov-Sahni braneworld cosmology 

The case A = — 1/(2|<t|) and v = 1 represents a class 
of braneworld inspired cosmologies due to Shtanov and 
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5 10 

log m^t 




5 10 

log m^t 



FIG. 1: Scaling solution for R-S II brane cosmology includ- 
ing a radiation fluid and a scalar field with potential given by 
Eq. 1501 . In order to obtain the scaling solution, we choose 
A = 10 and for simplicity specify the energy scale of the brane 
tension to be cr/m\ = 10~ 20 . In the upper figure, the time 
evolution of the energy density for both radiation p r (solid 
curve) and the scalar field (dashed curve) are shown. Note 
that around logm4i ~ 10 there is a rapid change as the 
quadratic correction becomes negligible and the standard cos- 
mological evolution is recovered. However, the energy density 
of the scalar field mimics that of the radiation fluid through- 
out the entire evolution, i.e., there is scaling behaviour. In 
the lower figure, we show the time dependence of the cos- 
mic expansion law. At early times, before logm4t ~ 10, the 
scale factor grows as a a t 1//4 , and represents a solution for 
a radiation-dominated universe in a p 2 dominated cosmol- 
ogy. After this time, the conventional expansion rate a oc t 1 ^ 2 
arises. 



Salmi (S-S) [25|, |39|, |4(J . In this scenario, a co-dimension 
one brane with negative tension a is embedded in a five- 
dimensional conformally flat Einstein space, where the 
signature of the fifth dimension is timelike. In this model, 
the deviation from the conventional Friedmann cosmol- 
ogy is characterised by 



i - 



(53) 



This model has recently been invoked to develop a non- 
singular oscillating universe, where the turning points in 
both the contracting and expanding phases are induced 
by the quadratic correction |42| . 

The S-S braneworld is dual to that of the RS-II sce- 
nario in the sense discussed above. The scaling potential 
follows directly from Eq. I|42|l : 



V(<p) = 2|cr|y; 2 sech 2 




(54) 



FIG. 2: As in Fig. Q, although now we specify A = 1. This 
leads to a scalar field dominated universe. In the lower figure, 
the time dependence of the cosmic expansion law is shown. 
As expected in a universe dominated by a scalar field, the 
scale factor grows as a oc t at early times, corresponding to 
the p 2 dominated phase, and as a oc t 2 at late times when the 
linear p term is important. 



whereas the time-dependences of the scalar field and 
scale factor follow from Eqs. (|48|l and l|47|) . respectively, 
after substituting for A = — l/(2|cr|) and v = 1: 



sinh 




Such a scaling solution is phcnomenologically interest- 
ing since it represents a non-singular bouncing cosmol- 
ogy. The universe collapses from infinity (t — » — oo) to a 
finite size at t = and then bounces into an expanding 
phase. The scalar field rolls up the potential during the 
collapse, reaches the maximum of the potential at cp = 
at the instant of the bounce, and then rolls monotonically 
down the other side during the expansion era. 



D. Cardassian cosmology 

In the above classes of models, the modifications to the 
Friedmann equation become significant at high energy 
scales (early times). On the other hand, recent CMB and 
large-scale structure observations indicate that the uni- 
verse is entering a stage of accelerated expansion at the 
present epoch and a number of phenomenological models 
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have been developed in an attempt to provide a geomet- 
rical interpretation of these observations. In Cardassian 
cosmology pfil l43j| . for example, the modification term in 
the Friedmann equation is given by Eq. (|4U|I with A > 
and v = n |26l l43| : 



L(p) = a/1 + Ap n 



(57) 



and the present-day cosmological acceleration can be ex- 
plained even when the energy density is comprised of only 
ordinary matter sources if n < —1/3. The characteristic 
feature of this model, therefore, is that the modification 
term becomes significant at late times. 

Although no scalar field is present in the scenario con- 
sidered in Refs. [SEl, it is instructive to show that 
the equivalent background cosmology can be obtained 
from a model comprised of a barotropic fluid and a 
self-interacting scalar field. We see immediately from 
Eq. Ij41|l that the corresponding potential which provides 
the fixed point attractor solution is given by 



V{(j>) = y c 2 A--cosech- 




(58) 



and in Fig. [3] we demonstrate numerically how the above 
potential leads to the expected attractor solution. Al- 
though such a homogeneous solution is indistinguishable 
from the purely perfect fluid background, it is possible 
that the presence of a scalar field may modify the evo- 
lution of perturbations and the clustering properties of 
matter. In principle, this could result in potentially ob- 
servational signatures and a scaling solution of this 
type provides a framework for quantitatively investigat- 
ing the evolution of perturbations in these models. 

Finally, before concluding this Section, we illustrate 
the duality transformation that relates the Cardassian 
cosmologies with the S-S braneworld. Denoting the for- 
mer with a subscript '+' and the latter by '— ', we may 
substitute the form for &+(0) in the Cardassian scenario, 
as given by Eq. (|4^|) . into Eq. $E\i to deduce that 



V-(4>) oc secli psz * 




(59) 



This reduces to the S-S scaling potential Eq. H42|l (with 
v = 1 and A = A_) when pfi^A 2 . — —n and A + = A_/n. 
Moreover, in this case, it can be verified that the dual 
correction function satisfying Eq. (|37|l reduces to Eq. 
(l4*S|l with v = 1 and A = A_. 

We now proceed in following section to determine the 
scaling solutions in a further braneworld scenario, where 
the corrections to the Friedmann equation become sig- 
nificant at late times. 



VI. DVALI-GABADADZE-PORRATI 
BRANEWORLD COSMOLOGY 

The Dvali-Gabadadze-Porrati (DGP) braneworld sce- 
nario pTl |45| corresponds to a 3-brane embedded in 
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FIG. 3: Scaling solution for a Cardassian-type cosmology in- 
cluding a matter fluid source and a scalar field with a po- 
tential given by Eq. 15811 . We specify n = —0.5 and choose 
A = 10 in order to obtain the scaling solution. For simplicity, 
we set the energy scale where modifications to the standard 
scenario become significant to be A/m\ = 10 -10 . In the up- 
per figure, the time evolution of the energy density of matter 
p m (solid curve) and that of the scalar field are shown. 
When log irnt ~ 10, the time-dependence of the energy den- 
sity changes as the late-time modification term becomes sig- 
nificant. In the lower figure, we show the time dependence 
of the cosmic expansion law. For log 77147; ~ 10, we see that 
a oc t 2 / 3 , corresponding to the conventional matter dominated 
universe, whereas for logm4t > 10, the correction term leads 
to an accelerating universe, a oc 



flat five-dimensional Minkowski spacetime, where a Ricci 
scalar term is included in the brane action. All energy- 
momentum is confined to the brane since the bulk is 
empty. The modified Friedmann equation for the DGP 
model is given by 



r Zvn\ 



(60) 



where r Q = 7714/(2771,5) and m 5 is the five-dimensional 
Planck scale. In the DGP model, gravity behaves as four- 
dimensional Einstein gravity at short scales, whereas it 
propagates into the bulk at large scales. This induces 
corrections to the standard Friedmann equation at low 
energies and the parameter ro determines the scale at 
which these corrections become important. There are 
two inequivalent ways of embedding the brane in the bulk 
and this is reflected in the different choices of sign in Eq. 
|(SDJ. In this section, we refer to these as the (+) and (— ) 
branches, respectively. In subsequent expressions, where 
different signs may be taken, the upper case corresponds 
to the (+) branch. 

It proves convenient to express the Friedmann equation 
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(|60|l in the equivalent form 



H = 



2r 



where 



B 



3ml 



(61) 



(62) 



For the (— ) branch, the late-time attractor is de Sitter 
(exponential) expansion for any decreasing energy den- 
sity j^. For the (+) branch, on the other hand, expand- 
ing Eq. I|61|) as a Taylor series to lowest-order implies 
that H w 87rrop/ (37714). Modulo a rescaling of the four- 
dimensional Planck mass, this corresponds formally to 
the high-energy limit (p ^> 2a) of the R-S II braneworld 
(149)1 . Consequently, the e arly -time analysis of the latter 
model performed in Ref. |3l| is directly applicable to the 
late-time behaviour of this branch of the DGP model. In 
particular, we may conclude immediately that the poten- 
tial driving the scaling solution in this limit is the inverse 
power-law potential V oc 4>~ 2 . 

A direct comparison between Eqs. Q and i)61|) implies 
that the Friedmann correction function is given by 



L = 



1 



(63) 



In order to derive the scaling solutions, we define a new 
variable, 8: 



p = — sinh 2 9. 
H B 



(64) 



Substituting Eq. ^0^1 into Eq. lj2U|) then implies that the 
solution represents a scaling solution if 



8ttA 



d6 



cosh# 



(65) 



and the integral 
of the potential: 



777,4 j cosh 6 =F 1 

1 ill may be evaluated to yield the form 



2n\ 



771,4 



Vl + Sp=Fl 



sinh \/ Bp. (66) 



The corresponding time dependence of the scaling so- 
lution can also be determined. In terms of the variable 
(|64|l . the Friedmann equation 1)611) simplifies to 



1 
H(9) = — sinh 2 - 
r 2 



for the (+) branch, and 
H{9) 



■ cosh 



ro 



(67) 



(68) 



for the (— ) branch. Recalling that <fi 2 /p = 2X 2 and 
Pcf, = Qtjjcp for the scaling solution, it follows that the 
scalar field equation l|21|) transforms to 



n 



3 A? 



de- 



cotanh 9 



(69) 



after substitution of Eq. Substituting Eqs. l|H7l) 

and 1)68)1 for the (+) and (— ) branches, respectively, and 
evaluating the integral 1)69)1 then implies that 



3X1 



t = cotanh 1 yl + Bp 



Vl + SpTl' 



(70) 



VII. SUMMARY 

In this paper we have brought together a number of 
recent approaches to cosmology which involve modifica- 
tions of the Friedmann equation. By introducing the gen- 
eral function L(p) as the way of parametrising the modi- 
fication, we have been able to establish the conditions un- 
der which the new system enters scaling solutions. Con- 
sidering the case where the energy density is comprised 
of a scalar field and background barotropic fluid, we have 
obtained the general relationship that would have to be 
satisfied between the evolving scalar field and L(p). In 
particular, we have obtained the corresponding potential 
V(4>) which leads to scaling solutions and, for a rather 
general class of functions of L(p), we have shown that 
there exist dual solutions which also exhibit similar scal- 
ing behaviour. This has allowed us to relate solutions 
which would otherwise appear quite distinct, including 
those involving collapsing and expanding cosmologies. 

Moreover, the duality can directly relate singular and 
non-singular cosmologies. To illustrate this property, 
consider a particular scaling solution a+ (0) that is singu- 
lar in the sense that the scale factor vanishes at a+(0) = 
0. (The value of the scalar field can be chosen to be <\> = 
without loss of generality). Suppose, however, that the 
logarithmic derivative of the scale factor with respect to 
the field is non-zero at this point, dlna+/d(f>\o 7^ 0, and 
furthermore, that the scale factor is a monotonic func- 
tion with a finite first derivative for all physical (non- 
zero) values of the field. These properties are satisfied, 
for example, in the R-S II and Cardassian models. 

The qualitative behaviour of the dual solution, a_(^>), 
is then determined from Eq. 1)3 ip. If we define a new pa- 
rameter e((j>) = [dlna/rf0] _1 , the value of the scale factor 
is simply given by the area under the curve e+ (</)) , where 
the field evolves from zero to some value <j>. (We are as- 
suming implicitly that <p > and p < again without 
loss of generality). However, due to the exponential na- 
ture of Eq. I|31|) the initial value of the dual scale factor is 
non-zero and this results in a non-singular background. 

On the other hand, the time reversal of the seed solu- 
tion a+ (4>) would result in a collapsing dual model where 
the limits in the integral (|31() are taken from <f> to zero. 
Consequently, the dual solution can be analytically con- 
tinued through (j> = into a contracting phase. In this 
sense, therefore, any singular (expanding) scaling solu- 
tion satisfying the above (very weak) conditions can gen- 
erate a non-singular bouncing cosmology, where the lat- 
ter is associated with a combination of the seed solution 
and its time reversal. For fixed values of {p, A±}, the col- 
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lapsing phase of the bouncing solution will be unstable 
if the expanding phase is stable, and vice-versa. In prin- 
ciple, however, different seed solutions may be employed 
to generate distinct and stable collapsing and expanding 
branches that can be smoothly joined at the bounce. 

This opens up the possibility that such dualities will 
allow us to relate singular cosmologies to non-singular 
bouncing cosmologies, a topic presently of considerable 
interest in cosmology. 

Finally, we have argued that the type of correction 
given by Eq. i|4U|) arises in a number of particle physics 
motivated models. Further examples arise in the limit 
where Ap v 3> 1. In particular, the case v = —1/3 cor- 
responds to the high-energy limit of the Gauss-Bonnet 
braneworld j^. In this model, the R-S II scenario is 
generalized to include a Gauss-Bonnet combination of 
curvature invariants in the five-dimensional bulk action. 
More generally, effective Friedmann equations of the form 
H 2 oc p v , where v is arbitrary |4l|, can arise in mod- 
els based on Hofava-Witten theory compactified on a 
Calabi-Yau three-fold [23] • Generalized scaling solutions 
driven by corrections of this form were recently investi- 
gated for a variety of scalar field models |23l |. 



Future directions involving the use of the duality prop- 
erties of these models would include an extension of our 
analysis to negative potentials, thereby allowing us to 
link these classes of solutions with those arising in the 
cyclic /ekpy rot ic scenario 0|. On the other hand, as 
we have seen, the more general corrections proposed in 
[27ll45l ] that arise due to modifications of gravity on large 
scales can lead to an explanation of the present cosmic ac- 
celeration without introducing dark energy (4^ ■ It would 
be interesting to investigate the impact that the duality 
transformations we have described have on such a sce- 
nario. 
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